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Under the excitation of strings, the wooden structure of string instruments is generally assumed to undergo linear 
vibrations. As an alternative to the direct measurement of the distortion rate at several vibration levels and frequen- 
cies, we characterise weak non-linearities by a signal-model approach based on cascade of Hammerstein models. 
In this approach, in a chain of two non-linear systems, two measurements are sufficient to estimate the non-linear 
contribution of the second (sub-)system which cannot be directly linearly driven, as a function of the exciting fre- 
quency. The experiment consists in exciting the instrument acoustically. The linear and non-linear contributions to 
the response of (a) the loudspeaker coupled to the room, (b) the instrument can be separated. Some methodological 
issues will be discussed. Findings pertaining to several instruments - one piano, two guitars, one violin - will be 
presented. 



1 Introduction: a linear behavior? 

Non-linear phenomena such as jump phenomenon, hys- 
teresis or internal resonance appear when the transverse vi- 
bration of a bi-dimensional structure exceeds amplitudes in 
the order of magnitude of its thickness Q . 

For string musical instrument (violin family, guitars, pi- 
anos...), the soundboard is generally assumed to undergo lin- 
ear vibrations: the transverse motion w remains in a smaller 
range than the board thickness. For example in the case of 
the piano, w remains in a smaller range, even when the pi- 
ano is played ff in the lower side of the keyboard. Askenfelt 
and Jansson fT\ report maximum values of the displacement 
at the bridge w^ax ~ 6 • 10"^ m in the frequency range 80- 
300 Hz (Fig.[T]). This maximum value is less than 10"^ times 
the board thickness (around 8 mm). It can therefore be as- 
sumed that, to a high level of approximation, the vibration of 
the soundboard is linear. 

The purpose of this article is thus to quantify experimen- 
tally such linear approximation. An original vibro-acoustical 
method is presented in Sec.[2]to isolate the soundboard non- 
linearity from that of the exciting device and to measure it. 
Soundboard intrinsic non-linearities of one upright piano, two 
guitars and one violin are then quantified using this method 
and results are presented in Sec. [3] 

2 A chain of two non-Unear systems 

When dealing with non-linearities, the non-linear contri- 
bution of the system under study (here a soundboard) has to 
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Figure 1 : Vibration levels at the bridge of a grand piano 
when played pp (dash-dotted line), mf (dotted line) and ff 
(solid line with • marks) for the notes C2 to B5 
(fundamental frequencies ^ 60 to 950 Hz), according to O. 



be quantified, and isolated from the non-linear contribution 
of the exciting device, which can be an electromechanical ex- 
citer or electrodynamic loudspeaker for example. The vibro- 
acoustic method presented in this section faces this problem 
and allows the estimation of the non-linear contribution of 
the soundboard of the instruments in cases where it cannot 
be directly linearly driven. 

2.1 Notations 

In Fig. [2] two non-linear systems modelling the exciting 
device (an electrodynamic loudspeaker) and the system un- 
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Figure 2: A chain of two non-linear systems. 



der study (the soundboard) are chained. The system Sp trans- 
forms its input signal X(/) into F(/), which becomes the 
input of the system Sg and Z(/) denotes the output of the 
whole chain. Sh stands for the non-linear system equivalent 
to the chaining ofSp and So- 
li is assumed here that only X{f) can be directly linearly 
driven and that both Y{f) and Z(/) are measurable. The con- 
tribution of the non-linear system Sq, which cannot be di- 
rectly linearly driven, has to be estimated. 

2.2 Modelling non-linear systems 

Volterra series are a convenient tool to express analyti- 
cally the relationship between the input eit) and the output 
^(0 of a weakly non-linear system 1 3 , 4 1 which is fully char- 
acterized by the knowledge of its Volterra Kernels in the fre- 
quency domain [Vkifi, ...Jk)}kew- 

Cascade of Hammerstein models constitute an interest- 
ing subclass of Volterra systems whose Kernels possess the 
following property: 



V/:,By,: V(/i,...,/,), y,(/i,...,/,) = Wi+...+/,) (1) 

Volterra Kernels of cascade of Hammerstein models can 
thus be expressed as functions of only one frequency variable 
and are in practice easy to estimate experimentally (H O. 
This simple method is based on a phase property of expo- 
nential sine sweeps and the Kernels of such a model can be 
estimated from only one measured response of the system. 

2.3 Non-linear system equivalent to the chain 

Let the Volterra Kernels {Fk{fi , fk)}kew , [Gkifi , fk)}kew 
and {Hkifi, fk)}kEW describe the systems Sf, Sg et, Sh- 
In the cascade-case presented in Fig. [2j Volterra Kernels of 
Sh can be expressed analytically as functions of the Volterra 
Kernels of Sf and Sg following ITj. For k = I, one obtains: 



Hi(fi) = F,{f,)Gi{fi) 



(2) 



This proves rigorously the intuitive result that the linear 
transfer function of a cascade of weakly non-linear systems 
is the product of the linear transfer functions of each systems 
composing the cascade. For k = 2, the following expression 
is obtained: 



Wl,/2)= F2(/l,/2)Gi(/i+/2) + ... 

... F,(f,)F,(f2)G2(fuf2) 



(3) 



Assuming that Sf and Sg can be modelled as cascade of 
Hammerstein models, Eq. ^ becomes: 



H2(fuf2) = F2(fl + /2)Gl(/l + /2) + . . . 

... ^l(/l)^l(/2)G2(/l+/2) 



A chain of cascade of Hammerstein models can thus in 
general not be modelled as a cascade of Hammerstein models 
as it does not fulfil Eq. ([T]) for ^ = 2 due to the ^i(/i)^i(/2) 
term. 

2.4 Non-linear contribution of the system Sg 

The output of the system Sg can generally be decom- 
posed in its linear and non-linear parts as follows: 



Z(f) 



(5) 



Gi(f)Y(f)^Z^^(f) 

Now suppose that Sf and Sh can be modelled as cascade 
of Hammerstein models and that their Kernels have been es- 
timated using the method proposed in (O. Linear transfer 
functions Fi (/) and Hi (/) and the signals Y(f) and Z(/) are 
thus known. The non-linear contribution of Sg is then: 



CSa(f) = 



NL 



if) 



(6) 



Using, Eqs. ^ and ([sjl and multiplying by Fi(J), Cscif) 
can then be conveniently computed as: 



CSaif) = 



Fi(f)Z(f) - mf)Y(f) 



Hi(f)Y(f) 

2.5 Numerical validation 



(7) 



The procedure described previously to estimate the non- 
linear contribution of the system Sg in the chain of Fig.[2]is 
now validated on a numerical example. 

Systems Sf and Sg have been modelled as cascade of 
Hammerstein models of order 4, and each of their Kernels 
as ARM A filters having two zeros and two poles. The poles 
and zeros of the diff'erent ARM A filters are given in Tab. [T] 
The amplitudes of the Kernels of Sf and Sg are presented in 
Fig. [3] The sampling frequency is chosen as equal to 96 kHz. 
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Table 1 : Poles and zeros of the ARMA filters used to 
simulate the non-linear systems Sf and Sg. 



The Kernels of the non-linear systems Sf and Sh have 
been afterwards estimated using the method presented in fSl 
between 20 Hz and 10 kHz, with 10 second exponential sweeps, 
and assuming non-linear systems of order 4. 

C^J^(f), the non-linear contribution of Sg is then esti- 
mated using Eq. ^ and its real value C^^^\f) is also com- 
puted from the knowledge of the diff'erent Kernels of the non- 
linear systems. C^Jl^(f) and C^^^\f) are plotted in Fig.|4j 

The agreement between the estimated non-linear contri- 
bution and the real non-linear contribution is found to be very 
good even if the system Sh is approximated by a cascade of 
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Figure 3: Amplitude of the Kernels of the simulated 
non-linear systems Sp (left) and of Sg (right). 



30 




Frequency (in Hz) 

Figure 4: C^J\f), the non-linear contribution of Sg 
estimated using Eq. (m) and its real value C^^^\f) computed 
from the differeniKernels of the non-linear systems. 

Hammerstein models which is mathematically not true here, 
as shown by Eq. Q. This thus validate the use of the pro- 
posed method to estimate the non-linear contribution of the 
system Sg which cannot be directly linearly driven in a chain 
of two non-linear systems. 

3 Application to string instruments 

3.1 Experimental protocol 

The soundboard non-linearities of four strings instruments 
of no particular merit have been estimated: one piano, two 
guitars and one violin. The experimental protocol (first de- 
veloped in 1 8 1), similar for each of the instruments, is pre- 
sented hereafter and drawn in Fig. [5] 

The instrument tuned normally and in playing conditions 
is put in a pseudo-anechoic room (anechoic walls and ceiling, 
ordinary ground). In the case of the guitars and the violin, 
the instrument is suspended vertically, clamped at the neck. 
A particular attention is taken to mute the strings by strips of 
foam (or woven in two or three places) inserted between them 
(see Fig. |6] for example). Two configurations [loudspeaker, 
room} and {loudspeaker, instrument, room} have been anal- 
ysed with the following procedure. The electrical excita- 
tion of the loudspeaker was an exponential swept- sine [50- 
4000] Hz (26 s duration). For each instrument, the amplitude 
of the loudspeaker was adjusted at the beginning of the study 




(a) (b) 



Figure 5: Experimental protocol: (a) first configuration 
(system Sp)- (b) second configuration (chain of systems 
Sh)' In the case of the piano, the vibrometer was replaced 
in the second configuration by accelerometers on the 
soundboard. 



so to obtain displacements of the soundboard correspond- 
ing to realistic playings. To give an idea, the G = 0.5 gain 
(see below), generates soundboards displacements at 500 Hz 
of around 10"^ m for the guitars and violin. For the pi- 
ano soundboard the displacement at 500 Hz is approximately 
10"^ m. According to Askenfelt and Jansson 1 2 1 these values 
correspond to the ff playing (see measurements at the bridge 
of a grand piano reported in Fig.[T]). 




Figure 6: Typical measurement on a violin (in playing 
conditions, with muted strings) excited by a loudspeaker. 
The reflective adhesive tape (where the velocity is measured 
by the laser vibrometer) is visible on the violin soundboard, 
just above the bridge. 

In the first configuration - {loudspeaker, room} - the acous- 
tic response of the room ^(0 is measured with a microphone 
placed in front of the loudspeaker (where the instrument is to 
be put in the second configuration). In the second configu- 
ration - {loudspeaker, instrument, room} - where the instru- 
ment replaces the microphone, the motion of the soundboard 
z(t) was measured with a laser vibrometer for the violin and 
the guitars (velocity in this case), and with accelerometers 
for the piano (acceleration in this case). 

Hence, the corresponding input-output scheme for the 
diff'erent configurations can be summarised by the diagram 
Fig. |2j where each system is weakly non-linear. The first 
configuration (a) in Fig. [5] corresponds to the non-linear sys- 
tem Sf and the second one (b) to the chain of two non-linear 
systems Sh- Hence, the method presented in Sec. |2] allows 



the estimation of the non-Hnear contribution of the sound- 
board of the instrument - second (sub-)system Sg - which 
cannot be directly linearly driven. 

For the present purpose, it will be assumed that the instru- 
ment response exhibits the same amount of non-linearity in 
the second configuration as it would if it was excited by the 
acoustical field that creates y(t) at the microphone. In other 
words, we consider that the instrument behaves, as far as 
non-linearities are concerned, like a slightly non-linear (and 
localised) microphone replacing the true one. 

3.2 Results 

An example of measurement showing the temporal wave- 
form and the spectrogram of the soundboard velocity mea- 
sured with the laser vibrometer is given in Fig. |7] (here for 
the first guitar). Modes of vibrations of the soundboard are 
clearly visible when instantaneous frequency of the excita- 
tion signal approaches the modal frequencies: increase of 
velocity (visible in the waveform) results in an increase of 
non-linearities. 




Time (s) 



Figure 7: Spectrogram and temporal waveform of the 
soundboard velocity of the first guitar during a typical 
measurements. 



For the four instruments measured, the spectra of the lin- 
ear parts of the responses and the relative non-linear contri- 
butions of the responses separated as explained in Sec. [2] are 
shown in Fig. [8] Fig. [9j Fig. [T0| and Fig. [TT] (respectively: 
piano, first guitar, second guitar and violin). For each fig- 
ures, the upper plots correspond to the Fourier Transform 
of the linear impulse reponse (typical "spectrum") and the 
lower plots correspond to the separated non-linear contribu- 
tions (with the same x- and scale to allow proper compari- 
son between the four instruments). The contributions Csp of 
the loudspeaker are given in dotted lines (first two figures), 
and the contributions Csc for the soundboard of the instru- 
ments in solid lines. Except for the piano case where only 
one loudspeaker's gain was tested, measurements have been 
done with four or five diflferent loudspeaker's amplitude. 

3.2.1 Piano 

The non-linearity content which can be attributed to the 
piano soundboard appears to be contained within -30 to -50 dB 
(see Fig. [8]). The apparent increase in non-linearity for fre- 
quencies below 100 Hz is probably an artefact of the method 
since the quality of the reconstruction of the non-linear im- 
pulse responses is degraded near the lower and upper bounds 
of the explored frequency range ([50-4000] Hz in the present 



case). Moreover, the increase of non-linearity at anti-resonances 
(see vertical green double arrows at 129 Hz and 219 Hz for 
example) is certainly due to the decrease of linear contribu- 
tion which enhances the measurement noise by construction 
(see Eqs. |6] and [7]). Altogether, the order of magnitude of 
-40 dB can be retained for the non-linear contribution of the 
piano soundboard at the level of playing. 
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Figure 8: Top - Fourier Transform of the linear impulse 

reponse of the piano soundboard (typical "spectrum"). 

Bottom - Measured relative non-linear contributions for 
both systems Sf and Sg- Dotted line: non-linearities of the 

loudspeaker C^^; Solid line: non-linearities of the piano 
soundboard Cs^ • (Linear contributions are equal to dB on 
this graph). 



Note that in the method proposed in this paper a high 
quality loudspeaker is essential for an efficient separation of 
the non-linearity contributions. This quality may be char- 
acterised by the amount of "acoustical" non-linearity in the 
first configuration (dotted line in Fig. [8]). For the piano mea- 
surements the non-linearities of the loudspeaker (Bose - 802 
Series II) are contained within -40 to -60 dB. In the guitars 
and violin cases (next three figures), the loudspeaker (Tapco 
S8) has non-linearities contained within -30 to -50 dB for 
the G = 0.5 gain. 

3.2.2 Guitars 



The results on the two guitars are given in Figs. [9] and 10 
for several gain G of the loudspeaker. The spectra of the lin- 
ear impulse responses are independent of the gain (almost 
undistinguishable one from the other; see the top of the fig- 
ures) which is consistent with the theory. 
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Figure 9: Non-linearities generated by the soundboard of 
the first guitar and estimated for different gain of the 
loudspeaker. Caption similar to Fig. [8] 
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Figure 10: Non-linearities generated by the soundboard of 
the second guitar and estimated for different gain of the 
loudspeaker. Caption similar to Fig. [8] (Dotted Hnes are 
removed for more clarity.) 
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Figure 1 1 : Non-linearities generated by the soundboard of 
the violin and estimated for different gain of the 
loudspeaker. Caption similar to Fig. [8] (Dotted lines are 
removed for more clarity.) 



As explained above, some peaks of non-linearity are clearly 
visible at anti-resonances, specially before the Helmholtz mode 
(first resonance) of the guitars (see green double arrows at 
71 Hz for the first instrument and 79 Hz for the second one). 

In the low and mid-frequency ranges the evolution with 
frequency of the non-linearities are similar for both guitars: 
a constant increase from -50 dB (at 100 Hz) to -30 dB at 
600 Hz. For frequencies higher, the non-linearities of the 
second guitar seems to reach a plateau level of -20 dB (for 
the two highest excitation level) whereas for the first guitar 
this increase continues and reaches -10 dB at 3.5 kHz. 

Moreover the non-linearities of the soundboards increase 
with G for both guitars, particularly in the frequency band 
[100-700] Hz which is also consistent with the theory. The 
fact that in the high frequencies this evolution of non-linearities 
with the gain is inverted for the second guitar may still be 
attributed to an increase of the signal to noise ratio (SNR) 
which leads to a decrease of the artefacts caused by noise 
(when G increases). This explains also why at anti-resonances 
the lower gains (lower SNR) increases the amount of esti- 
mated non-linearities (see for example the blue curves higher 
than black ones at anti-resonance of the Helmholtz modes, 
for both guitars). 

3.2.3 Violin 



Results on the violin are given in Fig. 1 1 Except for the 
peaks at anti-resonances the violin soundboard intrinsic non- 
linearity is contained within -20 dB to -50 dB. 

3.3 Discussion 

The assumption of linearity of the soundboard vibrations 
of the four instruments is verified to a high level of approxi- 



mation (see Fig. 12 where the results are plotted in the same 
graphs). The mean value of non-linearities is more than 20 dB 
less than the linear parts in all the cases, except near the lower 
and upper bounds of the explored frequency range ([50-4000] Hz) 
(artefact of the method) or at anti-resonances where the SNR 
decreases. The comparison reveals that the non-linearity con- 
tent of the first guitar (red line) is almost 10 dB more than for 
the other instruments, for frequencies higher than 1 kHz. A 
conclusion on the quality of this instrument is unfortunately 
impossible here (and is not the object of this paper). Such 
conclusions would require measurements with more instru- 
ments of different qualities/origins... 
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Figure 12: Comparison of the non-linearities generated by 
the four instruments. (G = 0.5 for the guitars and the violin.) 



4 Conclusion & Perspectives 

In this article an original vibro-acoustical method is pre- 
sented to isolate the soundboard non-linearity of string in- 
struments from that of the exciting device (here a loudspeaker) 
and to measure it. For a chain of two non-linear systems, 
the method allows the estimation of the non-linear contribu- 
tion of the second system which cannot be directly linearly 
driven. Experimental quantifications of the linear approxi- 
mation of the intrinsic soundboard vibrations of one upright 
piano, two guitars and one violin is given for level of exci- 
tation corresponding to thcj^ playing. These non-linearities 
are contained within 20 dB to 50 dB less than the linear parts 
in the [50-4000] Hz frequency range, except in the high fre- 
quency domain of one of the guitars. The measurement noise 
appears to be crucial for a proper estimation of these non- 
linearities. A technique which allows to separate the mea- 
surement noise from the system non-linearity and to measure 
it (as in [9]) must be integrated in the method before giving 
conclusions in terms of musical quality of the measured in- 
struments. 
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